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Abstract 



Exact analytical, closed-form solutions, expressed in terms of special functions, are presented 
for the case of a three-dimensional nonlinear quantum oscillator with a position dependent 
mass. This system is the generalization of the corresponding one-dimensional system, which 
has been the focus of recent attention. In contrast to other approaches, we are able to obtain 
solutions in terms of special functions, without a reliance upon a Rodrigues-type of formula. The 
^ ■ wave functions of the quantum oscillator have the familiar spherical harmonic solutions for the 

^1^ , angular part. For the s-states of the system, the radial equation accepts solutions that have been 

i~^^ ' recently found for the one-dimensional nonlinear quantum oscillator, given in terms of associated 

fT^ I Legendre functions, along with a constant shift in the energy eigenvalues. Radial solutions are 

l' ' obtained for all angular momentum states, along with the complete energy spectrum of the 

f— ^ , bound states. 
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1 Introduction 

An interesting one-dimensional model of a nonlinear quantum oscillator has been studied recently 
by Carihena, Rahada and Santander [4]. This model presents two interesting features, namely non- 
linearity of the potential V{x) and the presence of a position dependent mass M(x), producing a 
position dependent spring "constant" K{x) oc Af (x) in the oscillator potential. This feature allows 
the existence of solutions for the classical nonlinear oscillator that have a simple harmonic form. 
The quantum mechanical version of the model admits exact solutions for the wave functions and the 
energy eigenvalues, which were obtained by the authors of Ref. [3]. This was achieved through the 
use of a factorization method, and solutions for the wave functions could be represented by means 
of a Rodrigues-type of formula. The intriguing features and the solvabilty of this model makes it 
an attractive instructional example that could yield insights into other physical systems and math- 
ematical methods. Consequently, this same one-dimensional nonlinear quantum oscillator model 



was further studied by Schulze-Halberg and Morris |S], where the exact solutions are expressed 
in closed-form in terms of special functions, without the need to use the factorization method or 
Rodrigues-type of formula of Ref. [4J. The full solution spectrum is obtained, and is seen to exhibit 
orthogonality and normalizability, with a proper reduction to the harmonic oscillator model in the 
appropriate limit. 

Here we focus attention on the three-dimensional generalization of the ID model of Refs. [1],[9]. 
We note that this 3D version of the model takes into account not only vibrational degrees of 
freedom, but rotational ones, as well. That is, we consider all possible angular momentum states, 
not only s states which have a vanishing centrifugal term in the effective potential for the radial 
motion. The resulting energy spectrum therefore involves both a vibrational quantum number n 
and an angular momentum quantum number L. The effective potential for this 3D system is of a 
more complicated form, but has spherical symmetry, and presents a more realistic representation 
of a system with variable mass M{r) that can both vibrate and rotate. 

Position dependent masses can arise in condensed matter settings, as mentioned in Ref. [4J, and 
in theories of gravitation, such as scalar-tensor theories (see, for example, Ref. ^), where a fixed, 
constant mass in one conformal frame becomes a position dependent mass in a different conformal 
frame, as originally pointed out by Dicke [2]. An example of this, analyzed by Dicke [2], is seen in 
the prototypical Brans-Dicke theory [3j , where a test particle with constant mass ttiq in the "Jordan 
frame" representation of the theory has a position dependent mass m,[((/>(x^)] = A[(l){x^)]mo in the 
"Einstein frame" representation. Here, the function A{(f)) of the scalar field (j){x^) connects the 
metrics in the two different conformal frames: g^y(x") = A'^{(p)gfi^{x°'), with g^j^^ representing the 
Jordan frame metric and g^u the Einstein frame metric. Since the scalar field cp is, in general, 
a function of spacetime coordinates, the mass of the test particle in the Einstein frame becomes 
spacetime position dependent. Therefore, well motivated reasons exist for considering theories with 
position dependent mass. 

Here we consider the nonlinear quantum oscillator studied in Refs. [4] and [9], generalized to three 
dimensions. We follow the procedure of Ref. [4J of formulating the classical Hamiltonian in such a 
way as to allow us to pass to the quantum Hamiltonian, with a quantization procedure resembling 
canonical quantization. We arrive at the 3D time independent Schrodinger equation, which we 
separate into angular and radial parts. The angular solutions are the usual spherical harmonics. 
We then focus on the more difficult problem of solving the radial equation. This is achieved by 
presenting solutions in the form of special functions. The entire bound state energy spectrum is 
obtained in terms of vibrational and rotational quantum numbers. 

2 The classical system and its quantization 

In this section we derive a quantum model featuring a nonlinear oscillator potential in one and three 
dimensions. To this end, we start out from the classical Lagrangian and construct the quantum- 
mechanical Hamiltonian. In addition, we will briefly comment on solutions of the classical system. 



2.1 The one-dimensional case 

Our starting point of the analysis carried out in this note is the fohowing, one-dimensional 
oscillator potential that was first studied in [4]: 

where x is a real, scalar variable and m,a > stand for constants, m representing a mass parameter. 
The classical Lagrangian L for the system associated with the potential ([T]) is given by 
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note that the dot stands for differentiation with respect to time and that the abbreviation g = rac? 
was introduced. The Lagrangian (l2|) is seen to have a noncanonical kinetic term, and we identify a 
position-dependent mass M = M(x), related to the mass parameter m by 

^ m ^ g 

(Ax2 + 1) a2(Ax2 + l)- ^^ 

Let us note that the equation of motion following from our Lagrangian ([2]) is 

{\x^ + \)x- (A x) i^ + q2 X = 0. (4) 

Now, in order to obtain the classical Hamiltonian for the system governed by the Lagrangian ([2]), 
we extract its kinetic energy T: 

■1 1 2 

^ - M x^ = -!-—, (5) 
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where p denotes the canonical momentum, given by p = Mx. This leads us to define a modified 
momentum variable P, determined by the condition 

' ' (0) 



This relation renders the kinetic energy ([5]) in the form 

P2 p2 

T = -— = ^ (7) 

2m 2M' ^' 

such that the classical Hamiltonian for our system can now be written down explicitly as the sum 
of kinetic energy ([7]) and potential 

H = —P^ + V (8) 
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Note that the potential V is defined in dTJ. 



^ (A x^ + 1) p2 + V. (9) 



Quantization: The transition to tlie quantuni-meclianical context is done by introducing the 
momentum operator p = —i h d/dx, which converts ([5]) into the Hamiltonian operator H of our 
quantum system: 
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note that in the last step we have used the explicit representation of our potential ([T]). The 
Hamiltonian presented in Ref. 0] (see above Eq(6) of that paper) coincides with our finding (jlOp . 



2.2 The three-dimensional case 

We can extend the previous to three dimensions in a straightforward way, replacing the scalar 
variable x by a coordinate triple, denoted by r. We will work in spherical coordinates {r^6,ip), 
defined in the usual way, such that the classical Lagrangian ([2]) is converted its three-dimensional 
counterpart as follows: 
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where the abbreviation v = f was used and || • || stands for the euclidean norm in R^, note further 
that ||r|| = r. The position-dependent mass constructed in ([3]) maintains its form, except for the 
change of variable: 
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The equation of motion generated by the Lagrangian (jlip is very similar to its one-dimensional 
counterpart (JH): 
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Following the procedure shown for the one-dimensional context, we incorporate the mass condition 
([6]) and arrive at the following three-dimensional classical Hamiltonian 
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where V is taken from the exphcit form ([T]) with x^ replaced by r^. As a consequence, the potential 
is now spherically symmetric. 

Quantization: As before, the quantum-mechanical Hamiltonian operator H is obtained by ex- 
pressing the classical momentum P through the operator p = —i hV. This renders the Hamiltonian 
operator H in the following form: 

H = — P P + V 
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where the explicit form of V was taken from (|13p . 

2.3 Classical Solutions 

Before we start analyzing the three-dimensional quantum model, let us briefly comment on the 
solutions of the classical equations of motion. 

The one-dimensional case. Our equation of motion (j3]) admits the following solution [6] [3]: 

X = A sm{LO t + (j)), (15) 

where A, uj and the phase (p are real-valued constants, related through the constraint 

The solution psp has the form of a conventional harmonic oscillator. Note that we can put the 
potential into harmonic oscillator shape: from ([I|), ([3]), and ()16p we can rewrite the potential as 

V = - M a^ x^ 
2 

= ^ M{l + XA^) uj^ x^ 

= Ikx', (17) 



where we have employed the abbreviation 
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Note that in the last step we used relation ^, where x is replaced by r. Expression K could be 
seen as an effective spring "constant". 

The three-dimensional case. The equation of motion (jlip admits a solution that resembles 
its one-dimensional counterpart (|15|) : 

f = A sm{u} t + (/)), (19) 

where ^ is a constant amplitude, uj denotes the oscillation frequency and </> is a phase. As in the 
previous case, the latter parameters are subject to the constraint of (J16p . with A = \\A\\. This 
is not surprising, since this solution represents a one-dimensional oscillation with respect to some 
fixed axis in the {x, y, z) space. One can rotate the coordinate system so that the motion takes 
place along the x axis, with y = z = 0, in which case the three-dimensional solution degenerates 
to the single-variable situation. The function (J19p simply represents an s-state where the oscillator 
undergoes linear motion with no rotation. 

If rotational motion is taken into account, the situation is more complicated. Let us assume motion 
with vibration and rotation within the x-y plane, and use polar coordinates (r, 6) to describe the 
motion. We make use of the relations arising from twofold differentiation of r with respect to time 
t: 

f = r r v = rf + r96 a = {r — r 6 )f + {r 6 + 2 r 6) 9, (20) 

where r and 9 stand for unit vectors in polar coordinates, that is. 
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(21) 



Our equation of motion pip now yields that rO + 2r9 = 0, which implies r'^9 = C = constant. 
Using this, the radial part of equation (jlip can be written as 

(Ar2 + l)(r-^)-Ar(r2^^)+a2r = (22) 

This coincides with the one-dimensional equation of motion Q for zero angular momentum C = 0, 
that is, 9 = 0. However, for C ^ the motion is more complicated. 

3 The three-dimensional quantum system 

In the single- variable case, the quantum Hamiltonian (jlOp was shown to admit a discrete spectrum 
and an orthogonal set of corresponding eigenfunctions, both of which were obtained in closed form 
[H[9]. In this section we will show that the same is possible in the three-dimensional situation, 
where the dynamics is governed by the Hamiltonian p4p . 



3.1 Separation of variables 



We start out by setting up the stationary Schrodinger equation that is generated by our three- 
dimensional Hamiltonian (jl4[) . To this end, let us recall that the corresponding potential is ex- 
tracted from ([1]), where x^ is replaced by r^: 
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where g = m o? . This potential enters in our Schrodinger equation, which reads 
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Note that the solution ^ is allowed to depend on all three spatial variables. Now, the usual process 
of variable separation in spherical coordinates [7] renders \I' as a product of a radial and an angular 
solution, that is, ^(r, (/;,0) = R{r^ Y{lp,9). Taking the separation constant as L(L + 1) for a 
nonnegative integer L, the two functions R and Y are found to satisfy the following equations: 
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where the prime denotes differentiation with respect to the radial variable r. The two-dimensional 
angular equation (j26|) can be solved in terms of spherical harmonics [7|, usually denoted by y = 
Yj^(9,ip), where the integer /i satisfies | /x |< /. Since the spherical harmonics are well known and 
have been studied thoroughly, we will not comment on them any further. Instead, we focus on the 
radial equation (j25p and its solutions. 



3.2 The Radial Equation 

Before we start constructing the solutions of our radial equation ()25p . we will compare it to the 
one-dimensional Schrodinger equation studied in [H [9] , and afterwards employ a change of variable 
to rewrite it in dimensionless units. 



3.2.1 Comparison to the one-dimensional case 

In this paragraph we will show that the radial equation (j25p can be cast in a form similar to that of 
its one-dimensional counterpart. To this end, we need to gauge away the 2/r term in the coefficient 
of R' by introduction of a new function u = u{r), defined as 
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r R, 
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which obeys the fohowing equation 
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Now define a new energy parameter £ which is shifted from the energy E by an amount controhed 
through the parameter A: 



and define an effective potential 
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After substitution of the redefined energy ()29|) and potential ()30p into our equation ()28p , we obtain 
its simplified form as 
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More explicitly, using (p3]) 
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Comparing this to the one-dimensional Schrodinger equation of ^, we see that it is the same 
equation (with x replaced by r) for s-states, where L = 0. The corresponding solutions for the 
latter case, and as such for our radial equation in the form (|32p for L = 0, are given by associated 
Legendre functions [9] . Note also that the energy gets shifted by a constant as displayed in (j29p . Let 
us make one more remark on the effective potential (|30p. After incorporation of (|23p and g = ma^ 
we obtain 



V, 
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note that in the last step we used relation ^ with x replaced by r. We observe that expression 
(j33p exhibits an interesting symmetry, being invariant under the simultaneous transformations 
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analysis of which is beyond the scope of this note. 



3.2.2 A dimensionless radial equation 



We will now convert equation (j25p into a dimensionless form, which makes subsequent calculations 
much easier and more transparent. To this end, we define a dimensionless radial coordinate y and 
a dimensionless parameter A, such that Ar^ = Ay-^. The coordinates and parameters are related by 



^ = A 

C2 



Cy A = ;^ (35) 

where the constant C is to be determined. Using this, we rewrite Eq. (l25p in the form 
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Here it is understood that R is expressed through the new variable y, note that for reasons of 
brevity we have not changed the function's name. We now look at the term in ()36p that contains 
the potential F as a factor. After incorporation of (j23p . the latter term reads 



2mC2 
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Once we construct the solutions to our radial equation, it will prove useful to have a redefined 
constant g. Adopting the approach taken in [3], we set 
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where we used our transformation (j35p . This renders the term (j37p in the following form: 
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Next, we require Cm a /h = 1, that is, the free constant C is chosen to be 
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such that expression ([55]) obtains a simple, dimensionless form: 
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In the final step we plug (j39p and ()40p into our rewritten radial equation (j36p . which becomes 
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where the abbreviation e = E/(ha) was employed. We observe that the radial equation in its form 
TTl) is dimensionless. 



3.2.3 Solutions of the radial equation 



In this section we will show that our radial equation (j4ip admits a discrete spectrum and a corre- 
sponding set of orthogonal solutions. In order to simplify our calculations, we choose a system of 
units where h= 1, m = 1, and we set the speed of light c = 1. In this system of units, the dimension 
of distance is the reciprocal of the dimension of mass, both of which are now dimensionless. Let us 
now define the following intervals D\ C M: 
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(42) 



The case A = corresponds to the well-known harmonic oscillator model, as can be inferred from 
inspection of ()4ip , for details the reader may refer to [7] . Now, the boundary- value problem for our 
radial equation (I4ip . defined on D^, takes the following form 



iAy^ + l)R"+(- + 3Ay]R' + 
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equipped with A-dependent boundary conditions: 
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Note that we have expanded the coefficient of R in (j4ip . Relations (j43p . (j44p form a spectral 
problem with e being the spectral parameter. In order to be physically meaningful, solutions of 
equation (|43p must be located in the following weighted Hilbert space 
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The weight function /x that appears in the latter definition is given by 
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Note that the numerator of this weight function is due to the measure in spherical coordinates. Let 
us finally remark that the internal product of the space (jl5]) . that induces its norm, is defined as 
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where the asterisk stands for complex conjugation. Equation ()43p can be solved exactly, its general 
solution has the form 
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ci F + C2 G, 
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where the two functions F and G are given by the following expressions 

F = yL {Ky'' + 1)-^ 2Fi{a- P,a + P,L + ^,-KyA (48) 

G = y-^-i(Ay2 + i)2k2Fi('a-L-i + /3,a-L-^-/3,-L+i,-Ay2V (49) 

Here, 2F1 stands for the hypergeometric function, the arguments of which contain the following 
parameters 

L 1 1 

a = - + 50) 

2 2 2 A ^ ' 

. _ VA + l-2eA + A^(L^ + L + l) 

Since L is a nonnegative integer and since the hypergeometric function behaves like a constant 
around the origin, the solution (I49p is in general singular there and as such does not satisfy our 
first boundary condition in (j44p . Consequently, we must discard the solution (j49p by setting C2 = 
in ()47p . Focusing on the remaining solution (I48|] . we first recall that the hypergeometric function 
behaves asymptotically like an exponential function. Since this would lead to an unbounded solution 
and to a violation of our second boundary condition in (j44p . we will convert the hypergeometric 
function into a polynomial by requesting its first argument to equal a negative integer or zero. More 
precisely, we impose the condition a — (3 = —n, where n G NU {0}. Solving the latter condition for 
e leads to a discrete set (e^) C M of energies: 

LA 3 

en = -2 A n^ -2 L A n-2 A n \-2n + L + -. (52) 

We will now substitute this relation into our parameters (I50p , (I5ip and plug them into the solution 
()48p . For the sake of simplicity, we will identify the functions F and R, since they differ only by 
an irrelevant constant. Thus, we obtain the following set of solutions (Rn), n G N U {0}, 

Rn = y''{Ay^ + l)~^2Fi(-n,n + L + l-j,L + ^,-AyA. (53) 

Due to its first argument being a nonpositive integer, the hypergeometric function degenerates 
to a polynomial. In particular, the solution ()53p can be rewritten as the set of functions (Rn), 
n G N U {0}, given by 

Rn = y^(Ay2 + i)-2kp]^+i-i-|)(i + 2A/), (54) 

note that P stands for a Jacobi polynomial and that irrelevant factors have been discarded. We 
will now show that the solution set (j53p is orthogonal and that its elements are normalizable. To 
this end, we distinguish the two cases of negative and positive A. 

First case: A < 0. It will prove convenient to work with the representation (j54p of our solutions, 
which for negative A can be rewritten as follows 



Rn = /(-|A|y2 + i)^p„i^+^'W-^y(i_2|A|y2)_ (55) 
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In the first step we will show that these functions fulfill our boundary conditions (|44p and are 
normalizable. To this end, we observe that the solutions (j55p are continuous and bounded on Da, 
such that the first boundary condition verifies. Next, it is straightforward to see that the functions 



(j55p always admit a zero at r = y^1/|A|, because the power of the second factor on the right-hand 
side is always positive. In conclusion, all functions ([55]) comply with our boundary conditions (jM]). 
In order to investigate normalizability, we notice that the corresponding norm || • 11^2 is obtained 
by integrating the following expression over Da: 

|i?„|2/X = li^nP ^ =, (56) 

V-|A| 2/2 + 1 

where the weight function fj, was defined in ()46p and has been adjusted to negative A. Since 
expression (I56p is continuous on D/<^, we must verify boundedness at the endpoints of the latter 
interval. Inspection of (j55p and (j56p shows that |-R„p/i behaves at y = like 2/^^+2^ g^^id as such 



stays bounded there. At the remaining endpoint y = y^l/IAI of our domain we find the asymptotics 

(— |A| y2 _|_ i^TO 2. We see that although the integrand (j56p becomes unbounded for A < —2, it 
remains integrable, because its exponent stays above -1. Next, in order to show that the solution 
family (Rn) is orthogonal with respect to the inner product in L^ (Da)) we make use of the following, 
well-known identity [l] 

1 

J{l-xr{l + xf Pt''Hx)Pt'\x)dx = C5mn, (57) 

-1 

valid for a, 6 G M, a,b > —1, and m,n G NU {0}. Furthermore, 5 denotes the Kronecker symbol 
and c G C is a constant, the explicit form of which is not of importance here. Let us apply the 
following settings to relation (|57p : 



a = L + ^ 6=i^-i x = l-2\A\y^, (58) 

It is immediate to see that a,b > —1, such that (j57p is applicable. Absorbing all constant factors 
into a parameter k £ C, the left-hand side of (|57p obtains the following form 

1 

-1 
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Combination of this result with (j^Tj) gives the identity 

{^Rmi Hn) [2 = C Omni 

where C = c/k G C. In summary, for A < 0, the solution of our boundary- value problem (|43p . 
(HH) is given by the infinite number of functions Rn € L^(L>a)> n € N U {0}, as displayed in 
(j^5|) . In particular, the set {Rn) is orthogonal with respect to the inner product of L'i[D\). The 
corresponding spectral values (e„), n G NU [{0} can be found in (j52]) : 

en = 2 |A| n^ + 2 L |A| n + 2 |A| n + — Li + 2 n + L+-, (59) 

where the substitution A = — |A| was made. 

Second case: A > 0. As in the previous case, we start out by verifying that our solutions ()54p 
satisfy the boundary conditions (j44p for the present case of positive A. Inspection of the explicit 
form (f54l) shows that the functions R^ n G N U {0}, behave like y^ close to y = 0, such that 
our first boundary condition in (j44p is fulfilled. In a similar manner we can find our solutions' 
asymptotics at infinity, which turns out to be 7-^-i/^+2"-_ jf w^is expression tends to zero, then 
our second boundary condition in (I44p is fulfilled, such that we are left with a constraint on our 
solution index n, that is, 

n < ^-| (60) 

This constraint already indicates that the number of physically meaningful solutions cannot be 
infinitely large, as in the previous case of negative A. Before we comment on this matter in more 
detail, let us verify normalizability of our solutions, assuming the above condition on n to be 
fulfilled. The norm || • ||l2 of our solutions is obtained by integrating the expression 

|i?„|V = |i?n|'^===, (61) 

V A y2 + 1 

over I?A) which in the present case is the positive real axis. Inspection of our solutions ()54p shows 
that (j6ip is continuous on Da, such that only the asymptotics of (j6ip at infinity can affect normal- 
izability of the functions (j54p . It is straightforward to verify that (|6ip behaves like y2L+i-2/A+4n 
at infinity, which results in the condition 

1 1 L 

n < — . 62 

2 A 2 2 ^ ' 



We see that this condition on n is stronger than its previous counterpart (I60p . such that we can 
work solely with (I62p . If the latter condition is satisfied, then the corresponding function (I6ip can 
be integrated over D\. In other words, Rn G L'i{D\) for all n that ([62]) applies to. Hence, in 
contrast to the previous case of negative A, the number of physically acceptable solutions to our 
boundary- value problem ([I3|) . (jM]) is limited and depends on A. In particular, if the relation 

A > ^, 

- 1 + L' 
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holds, then there is no physical solution at all. Now that we have completed our analysis on 
normalizability, it remains to show that the functions (j54|) form an orthogonal set with respect to 
the inner product on L^(Da)- To this end, we will work with the hypergeometric representation 
()53p of our solutions. Taking into account that the latter functions are real-valued, we obtain for 
n,m £ NU {0} and m ^ n, the internal product of two different solutions as 

{Rm,Rn)i2 = / Rm Rn fJ- dy 

Da 

oo 



X ^Fi(-n,n + L + l-j,L + ^,-AyAdy. (63) 

Note that our weight function (j46p has been absorbed into the first two factors of the integral. In 
the next step we will employ the following change of variable: 



1 — X , , 

It is straightforward to see that this function maps x G (—1, 1) onto y G (0,oo). The differential of 
is given by 



such that after substitution of the latter relation and (j64p . our inner product (j63p takes the following 
form: 

1 

{Rm,Rn)Ll = ai /"(l-x)^+5 (l + a:)-2-i+i 2i^i f-m,m + L + l-|-,L + ^,|^j X 

-1 

X 2Fi(-n,n + L + l--,L + -,^——-]dx. (65) 

Note that all irrelevant constant factors have been absorbed into ai G M. In the next step we 
convert the hypergeometric functions into Jacobi polynomials, using the following identity [8] 

,Fj-„.-n-i,a+u'-^) = (4,^)"i>±il!ip(-)M. (66) 

\ x—lj \x + lj r(n + a + 1) 

valid for n G N U {0} and a,b,x G M, provided the hypergeometric function exists. In order to use 
relation (I66p in the present case, we must choose a = L + 1/2 and h = —2n — L — 1 + 1/A resp. 
b = —2m — L — 1 + 1/A, as inspection of (|65p shows. The inner product then renders in the form 

1 

{Rm,Rn)Ll = a2 |(l-x)^+^(l + x)-2-^+*— "pi^+^'-'-^-""+^)(x)x 



-1 



X pi'^'^'-'-'-'"'^^^\x) dx, (67) 
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where constant factors from (|66|) have been put into 02 G M. We observe that the internal product 
()67p is invariant under exchanging m and n. Without restriction, we can therefore assume that 
m > n and replace the Jacobi polynomial of order m by its representation through Rodrigues' 
relation [1]: 



p{a,i 



^ ^^ {l-x)-'^{l + x)-' "^ 



(1 - xY+'^ (1 + xf^"" 



where we must use the parameter settings a = L + 1/2 and h = —2m — L — 1 + 1/A: 

1 



[Km, Hn) ip. 



0-3 



/' 



1 + x) 



-l+m—n 



(F' 



dx^ 



{i-xY 



X pi"+^'-^-"-^"+^)(x)d., 



(68) 



(69) 

note that constant factors from (|68|) have been placed into 03 G M. In the next step, we perform 
an n-fold integration by parts, which yields 



\^mi ^nj 



Ll 



«3 



[l-X 





"(l + ,)-l+r„-n^(^+i-l~^-2n+^)^^^- 


^L+i^ 


'"^ (l + x)-™-^-^+^ dx. 



(70) 



Let us point out that the intermediate terms of each integration by parts vanish due to nonnegative 
exponents. In particular, the last exponent under the derivative in ()69p is positive, because 



1 1 1 

m-L-l + - > -2m-L-- + - > 0, 
A 2 A 



due to our assumption ()62p . Our final step consists in observing that in (|70p the m-fold derivative 
yields zero, because it is applied to a polynomial of degree —1 + m — n + n = m — 1. Therefore, we 
conclude 

{Rm,Rn)L2 = as Omm 

which establishes orthogonality. In summary, for A € (0,oo), the solution of our boundary- value 
problem (03]), (03]) is given by the functions i?„ G L^(Da), n G NU {0}, as displayed in ([55]). The 
number of functions is limited and depends on A through relation (I62p . The set (Rn) is orthogonal 
with respect to the inner product of L^ (Da) and the corresponding spectral values (e„), n G NU[{0}, 
are given by (f52]) . 



3.2.4 Convergence to the harmonic oscillator 

We will now briefly verify that the solutions and spectral values found in the previous section reduce 
correctly to their well-known counterparts if our parameter A goes to zero. Since the angular part 
of our model does not involve A, it is sufficient to study the radial boundary- value problem (03])) 
(j44p and its domain D/i^, as given in (j42p . Starting out with the latter domain, we observe that 
its two distinct definitions for negative and positive A become the same in the limiting case, that 
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is, liniA-s^o-DA = (0,00). Applying the same limit to the radial equation (jl5|) . our boundary- value 
problem reads 



R" + - R' + 

y 



y 


R = 


= 0, 


R{0) < 00, lim R{y) = 

y^oo 


= 0. 



ye (0,00) (71) 

(72) 



This coincides with the well-known radial equation for the three-dimensional harmonic oscillator 
^. After having determined convergence of the boundary- value problem, let us now investigate its 
solutions. As far as the corresponding spectral values are concerned, their explicit forms (|52p and 
(j59|) tend to the same limit if A goes to zero: 

en = 2n + L + | (73) 

which gives precisely the energy spectrum of the harmonic oscillator in three dimensions [7]. In the 
final step we must apply the limit A — t- to our solutions, where we will take the representation 
([53|) . Using standard identities [1], we find 

lim y^ {A y^ + 1)-^ 2F1 ( -n,n + L + 1 - ^,L + ^,-A y^ 

= y^exp(-^yAiFJ-n,L+^,y'^ 

= y^exp(-iy2^ L^n^Hy'), (74) 

where some irrelevant factors have been omitted. Note that i-Fi stands for the confluent hyperge- 
ometric function and L represents a generalized Laguerre polynomial. The functions ()74p are the 
well-known solutions of the boundary- value problem (I7ip . ()72p that lie in L^(0,oo) and belong to 
the discrete spectral values given by ([73]) . 

4 Concluding remarks 

In this note we have studied the quantum model of an isotropic nonlinear oscillator in three spatial 
dimensions. Via separation of variables in the underlying Schrodinger equation, we obtained the 
entire discrete spectrum of the problem, together with a corresponding set of orthogonal solutions. 
Connections to the one-dimensional situation and the conventional harmonic oscillator system were 
established. A natural continuation of this work could concern the behaviour of our model in a 
setting of spherical symmetry and arbitrary dimension (hyperspherical symmetry). In such a case, 
properties of our system like the energy spectrum and the number of physically acceptable solutions 
become functions of the dimension. 
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